Four-atom period in the conductance of monatomic Al wires 
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We present first principles calculations based on density functional theory for the conductance of 
monatomic Al wires between Al(lll) electrodes. In contrast to the even-odd oscillations observed in 
other metallic wires, the conductance of the Al wires is found to oscillate with a period of 4 atoms as 
the length of the wire is varied. Although local charge neutrality can account for the observed period 
it leads to an incorrect phase. We explain the conductance behavior using a resonant transport model 
based on the electronic structure of the infinite wire. 

PACS numbers: 73.40.Jn, 73.40.Gk, 73.63.Nm, 85.65.+h 



Metallic chains of single atoms represent the ulti- 
mate limit of the miniaturization of electrical conductors. 
Apart from possible technological applications as inter- 
connects in molecular electronics, the simple structure of 
these monatomic wires makes them an ideal test ground 
for developing and validating our understanding of elec- 
tron transport on the nanometer length scale. The con- 
ductance of monatomic wires has been investigated the- 
oretically by several authors using a number of different 
techniques HEHSHH Lang was the first to study the 
dependence of the conductance on the number of atoms 
in the wire. In 1997 he found that the conductance of a 
chain of Na atoms between jellium electrodes oscillates 
with a period of 2 atoms as the length of the chain is 
varied 1], and the following year Lang and Avouris ob- 
served the same behavior for chains of C atoms 0] . This 
even-odd effect was also found in Na wires by Tsukamoto 
and Hirose [3j using a method similar to that of Lang, 
and by Sim et al. [j] who combined cluster density func- 
tional theory (DFT) calculations with the Friedel sum 
rule to obtain the conductance. Recently, the even-odd 
oscillations have been confirmed experimentally for Au, 
Pt and Ir by means of the mechanically controlled break 
junction 0. 

The fact that a large number of metals of rather dis- 
tinct character exhibit the even-odd effect has lead to 
the suggestion that this is a universal feature of atomic 
wires [7J. In this paper we present calculations based 
on DFT showing that the conductance of monatomic Al 
wires varies in an oscillatory manner as the number of 
atoms in the wire is changed. In contrast to the even-odd 
effect, however, the oscillation has a period of 4 atoms. 
Based on a resonant transport picture we show that local 
charge neutrality can account for the observed 4-atom 
period. For long wires, however, we find that the oc- 
cupancy of molecular resonant states and thus the net 
charge located on the wire is determined by the position 
of the discrete energy spectrum of the free wire relative 
to the Fermi level of the electrodes. In general these two 
effects compete in determining the phase of the conduc- 



tance oscillation. In the case of Al we find that the latter 
effect dominates for wires containing more than about 4 
atoms, whereas charge neutrality is more important for 
shorter wires. By relating the electronic structure of the 
free wire to that of the infinite wire, we arrive at a sim- 
ple picture for the position of the free-wire valence reso- 
nances which accounts for the conductance both qualita- 
tively and quantitatively. In this picture the filling factor 
of the infinite wire determines the period of the conduc- 
tance oscillation: a half filled band implies the even-odd 
effect while a filling factor of 0.25 leads to a 4-atom pe- 
riod. 

We represent the combined electrode-wire system in a 
supercell containing the wire itself together with seven 
3x3 sections of the (111) planes in bulk aluminum. An 
example of such a supercell is shown in Fig. ^ Peri- 
odic boundary conditions are imposed on the supercell 
in the directions perpendicular to the wire axis, and for 
the DFT calculations also in the parallel direction. The 
calculational procedure consists of two main steps: (i) 
first a realistic wire configuration is found by relaxing the 
ion positions as well as the electrode-electrode distance 
as determined by the length of the supercell. During the 
relaxation the wire is constrained to be linear and ions 
below the surface layer are fixed in the bulk crystal struc- 
ture. These simulations were performed within DFT us- 
ing a plane-wave based pseudopotential code @| with an 
energy cutoff of 15 Ry for the plane- wave expansion. The 
ion cores are replaced by ultrasoft pseudopotentials @ 
and to treat exchange and correlation we use the PW91 
functional The bond lengths in the wires lie in the 
range 2.36 A to 2.50 A with the shorter bond lengths 
achieved closer to the center of the wire. For the longer 
wires the bond lengths vary less and approach a value of 
2.39 A which is also found for the infinite wire. No dimer- 
ization is found in accordance wi th p revious results for 
infinite linear chains of Al atoms [ll|. We mention that 
Ayuela et al. have found that an infinite linear Al chain 
has a nonmagnetic groundstate for interatomic distances 
smaller than 2.7 A 12] . (ii) The self-consistent effective 
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FIG. 1: The N = 3 wire connecting Al(lll) electrodes via 
a 3-atom basis. This structure constitutes the supercell in 
the DFT calculations and forms the scattering region in the 
conductance calculations. 



potential generated by the DFT code is used as input to 
the Green's function transport code described in detail 
in Ref. 01 ■ The effective potential of the supercell forms 
the scattering region in the transport problem while the 
potential in the semi-infinite leads is obtained from a 
DFT bulk calculation. In this way the scattering region 
and the leads are treated on an equal footing taking their 
full atomic and electronic structure into account. The 
conductance is calculated from the Green's function of 
the scattering region which is represented in terms of a 
system-independent basis set containing wavelets of com- 
pact support. All conductance calculations have been 
converged with respect to the size of the basis set. Due to 
the limited size of the supercell in the directions perpen- 
dicular to the wire, the two-dimensional Brillouin zone is 
sampled by 8 special k-points, both in the transport and 
the DFT calculations. 

In Fig. |21 we have plotted the calculated conductance 
of Al wires containing 1 to 9 atoms and attached to the 
(111) surface cither directly or via a 3-atom basis. We 
regard an atom as part of the wire if it has a coordination 
number of 2, see Fig. ^ With this convention the wire 
couples to the electrodes via a single atom which we refer 
to as the contact atom. The conductance oscillates with 
a period of 4 atoms taking values in the range 0.5Go 
to 1.7 Go- The oscillation amplitude increases approxi- 
mately 10% when the 3-atom basis is used as contact to 
the electrodes, but otherwise no difference is found be- 
tween the two geometries. 

We can gain some insight into the electronic struc- 
ture of the current carrying states by studying the eigen- 
channel resolved conductance [1J| as shown in Fig. [3| 
for N = 3. The valence configuration of atomic Al is 
3s 2 3p, which for the infinite wire leads to a fully occu- 
pied cr-band and a degenerate 7r-band with filling factor 
/ = 0.25. Three channels contribute to the total con- 
ductance: two similar 7r-channels and a single cr-channcl. 
The two 7r-channels constitute the main part of the con- 
ductance, while the c-channel is a tunneling channel 
whose contribution decays exponentially with the wire 
length from - 0.2G at N = 1 to ~ 10~ 4 Go at N = 9. 

To understand the origin of the conductance oscilla- 
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FIG. 2: Calculated conductance as a function of wire length 
with and without a 3-atom basis as contact to the electrodes. 
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FIG. 3: Eigen-channel resolved conductance for the N = 3 
wire shown in Fig. Q The two similar 7r-channels together 
with the er-channel are indicated by the dashed lines, while 
the solid line shows the total conductance. 



tions we first relate the conductance to the electronic 
structure of the wire. When a free iV-atom wire is con- 
nected to bulk electrodes, the molecular orbitals on the 
wire hybridize with extended states in the electrode. As a 
result of this mixing, the n'th energy level of the free wire 
is broadened into a resonance described by the projected 
density of states 



Pn{e) = 



l A n (e) 



A n ( £ ))2+A£( £ ) 



(1) 



For each molecular level, the coupling to the electrodes 
defines a unique orbital, the group orbital, which is lo- 
cated on the contact atoms. In terms of the group or- 
bital we can write A n (e) = Tr\V n \ 2 p c ,n(£) , where V n is the 
matrix element coupling the n'th molecular level with its 
corresponding group orbital and p c >n is the projected den- 
sity of states for the group orbital in the absence of the 
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wire. For simplicity we shall not distinguish between the 
group orbitals and thus we set p c>n = p c . If we further- 
more neglect interference between the different molecular 
levels we obtain a simple expression for the conductance 

G=^]>> 2 |K| 2 Pc(eFK(£F), (2) 

n 

where the sum extends over all molecular levels. From 
this expression we expect the position of the resonances 
relative to the Fermi level of the electrodes to be a crucial 
parameter for the conductance. 

So far, local charge neutrality has been considered the 
main cause of the even-odd oscillations in monovalent 
wires 0, IE • Indeed, for monovalent wires consisting 
of an odd number of atoms, local charge neutrality pins 
the Fermi level to the middle of a resonance resulting in a 
high conductance. On the other hand, for even-numbered 
wires the Fermi level lies off resonance leading to a lower 
conductance. In the case of Al the degeneracy of the tt- 
orbitals on the free wire and the fact that each Al atom 
provides one electron to the 7r-system, implies that local 
charge neutrality would cause oscillations in the conduc- 
tance with a period of 4 atoms. However, the resulting 
phase of the oscillation would be incorrect. The maxima 
(minima) occur for half-filled (filled) resonances which 
should correspond to N = 2, 6, . . . (N = 4, 8, . . .)• These 
numbers are off by 1 atom as compared to Fig. [21 Con- 
sequently, charge neutrality alone cannot account for the 
conductance of the Al wires. At this point it is interesting 
to notice that different results have been obtained for the 
phase of the conductance oscillation in Na wires. As op- 
posed to Tsukamoto and Hirose and Sim et al., Lang jlj 
has found a higher conductance for even-numbered than 
for odd- numbered Na chains. Calculations by Havu et 
al. || for Na wires between jellium cones show that it 
is possible to change the phase of the conductance os- 
cillation from even-odd to odd-even by varying the lead 
opening angle. 

To explain the behavior of the conductance, we need 
to examine more closely the position of the resonances 
relative to the Fermi level of the electrode-wire-electrode 
system. As a first step in this direction we consider a lin- 
ear chain of N sites coupled by nearest neighbor hopping. 
The Hamiltonian is 

N N—l 

H = Y, £ ° c ^ + J2 M+x c i + c h+i) (3) 

i=l i=l 

where a annihilates an electron at site i. The eigenvalues, 
, of H can be expressed in terms of the band structure 
of the infinite chain of coupled sites, e(k), as follows: 

4 = e(^ I ),n = l,...,iV (4) 

where we have set the inter-site distance to 1. Although 
the above relation is only strictly valid for the model 



Hamiltonian @, we can still use it to estimate the spec- 
trum of a free iV-atom wire in terms of the band structure 
of the infinite wire. As already described, the coupling of 
the free wire with the electrodes converts the discrete en- 
ergy spectrum into a set of resonances. Assuming p c (e) 
to be constant the n'th resonance of the iV-atom wire be- 
comes a Lorentzian of width Ajv .„ = w\ Vjv.nl 2 Pc- Since 
we are only interested in resonances close to the Fermi 
energy, it is reasonable to assume an energy independent 
coupling, Vjv.n = Vn ■ Furthermore, the extended nature 
of the eigenstates of the free wire causes \Vn\ 2 oc 1/N. In 
Fig. 0] we have illustrated in /c-space the evolution of the 
discrete spectrum of a free wire as a function of the wire 
length. For each eigenvalue we have indicated the width 
of the corresponding resonance. The eigenvalues are con- 
verted to energy via the valence band of the infinite wire 
and relation (J3J. To find the position of the Fermi level 
in the electrode-wire-electrode system we note that the 
Fermi energy is fixed by the macroscopic electrodes, and 
thus independent of the length of the wire. Since we do 
not expect a wire to charge up as it grows longer the net 
charge per atom must vanish in the limit of a long wire, 
and the band filling, /, must approach that of the infi- 
nite wire. In Fig. 0] we have indicated the Fermi level for 
the Al wire by a vertical line at k = 7r/4 corresponding 
to / = 0.25. As N increases the resonances move down 
through the Fermi level in a systematic fashion. The co- 
incidence of the Fermi level with the center of a resonance 
is a periodic event of period 4 starting at N = 3. Like the 
local charge neutrality, this picture explains the observed 
period of the conductance oscillation but in addition it 
provides information about the actual position of the res- 
onances. The first maximum (minimum) is predicted to 
occur at N = 3 (N = 5), which indeed gives the correct 
phase as compared to Fig.|3 The vertical line at k = tt/2 
indicates the position of the Fermi level for metals with 
a half-filled valence band such as Na, Au and C. Clearly, 
our model accounts for the even-odd behavior observed 
in these wires. 

According to the filling of 7r-resonances as predicted by 
Fig. |H the AT-atom Al wire should hold approximately 
N—l valence electrons. For short wires the electro- 
static potential from the missing electron drags the reso- 
nances down in energy thereby increasing the occupancy. 
The resulting position of the resonances must be deter- 
mined self-consistently and will in general lead to an oc- 
cupation somewhere between N—l and N. For short 
wires the simple picture is therefore modified by a charge- 
neutralization process which tend to shift down the phase 
of the conductance oscillation by one atom. This ef- 
fect explains the relatively high conductances found for 
N = 1, 2 as compared to N — 5, 6. 

To test our simple model quantitatively we have used 
it to reproduce the first principles results. We use Eq. 
(J2J) for the conductance and assume a Lorentzian shape 
of the projected density of states. The molecular lev- 
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FIG. 4: Evolution in fc-space of the discrete spectrum of a free 
iV-atom wire. The interatomic distance in the wire is set to 1. 
The horizontal lines symbolize the width of the corresponding 
resonance formed when the free wire is coupled to electrodes. 
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FIG. 5: The conductance of Al wires as calculated from first 
principles (solid) and a simple model (dashed). 



els given in Eq. (J1J are obtained via the 7r-band of the 
infinite Al wire calculated within DFT. From the first 
resonance in Fig. we read off Ajv = 3 = 0.15 eV as the 
half- width at half-maximum. Due to the 1 /N scaling, all 
resonance widths then follow from Ajv=3- The result of 
the model is seen in Fig. [5] together with the first prin- 
ciples calculation, repeated here for convenience. Except 
for the shortest wires (N = 1,2) where charge transfer as 
well as conduction through the cr-channel is important, 
the model gives a good description of the conductance 
with an average deviation of less than 15%. The fact 
that the model gives values above 2Go is a consequence 
of the neglected interference between different molecular 
levels. Finally, we note that the l/7V-dependence of the 



resonance width together with the A^-dependence of the 
resonance spacing implies the existence of conductance 
oscillations also for longer wires. This is in agreement 
with previous results in Ref. [lfij . 

In summary, our calculations show that in contrast 
to the even-odd effect the conductance of monatomic Al 
wires oscillates with a 4-atom period as the number of 
atoms in the wire is varied. This behavior is explained 
by combining a resonant transport picture with the elec- 
tronic structure of the free Al wire. By relating the elec- 
tronic structure of a free wire to the valence band of the 
infinite wire, we find that the period of oscillation is de- 
termined by the filling factor of the valence band of the 
infinite wire. For long wires (> 4 atoms) a simple model 
based on these ideas was shown to reproduce the first 
principles conductances while for short wires the simple 
picture is modified by charge transfer and tunneling ef- 
fects. 
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